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Abstract. We present the results of a radio monitoring campaign to search for radio emission from nearby ultra-luminous 
X-ray sources (ULXs). These sources are bright off-nuclear X-ray point sources with luminosities exceeding L x > 10 39 erg/sec. 
A well-defined sample of the 9 nearest ULXs has been monitored eight times during 5 months with the Very Large Array in A 
and B configuration. Our limiting sensitivity is a 0. 15 mjy (4 <x) for radio flares and » 60 /jjy for continuous emission. In M82 
two ULXs seem to have coincident compact radio sources, which are probably supernova remnants. No continuous or flaring 
radio emission has been detected from any other ULX. Thus, ULXs do not generally emit steady-state radio emission above 
radio powers of 1.5 x 10 17 W/Hz. The non-detections of the continuous emission are consistent with beamed or unbeamed 
radio emission from accreting black holes of < 10 3 M o based on the radio/X-ray correlation. Other published radio detections 
(M82, NGC 5408) are also discussed in this context. Both detections are significantly above our detection limit. If ULXs have 
flaring radio emission above 4 x 10 17 W/Hz we can give an upper limit on the duty circle of the flares of 6%. This upper limit 
is in agreement with the observed number of flares in Galactic radio transients. Additionally we present a yet unreported radio 
double structure in the nearby low-luminosity AGN NGC 4736. 



1. Introduction 

Ultra-luminous X-ray sources are among the most widely dis- 
cussed objects, but their true nature is yet unknown. These ob- 
jects have X-ray luminosities Lx ~ 10 39 ~ 42 erg/s which seems 
too bright for normal (stellar-mass) black hole X-ray binaries 
(BHXRBs), but far dimmer than normal active galactic nu- 
clei (AGN). The first hints to these intermediate-lumi nosity 
X-ray point sources have been found in the 1980s ( Fab biand 
Il989t IColbert et all 1 19951) . Observations with ROSAT, and 
subsequent high resolution X-ray satellites, confirmed these 
findings and showed that these intri guing sources are of- 
ten not in the centers of the gal axies ( Colbert & Mushotzkv 
1999; Roberts & Warwick 2000). More recent surveys (see 
e.g.. lColbert & Ptakll2002l IPtak & Colbertll2004 find approx- 
imately one ULX in every five galaxies, confirming that ULXs 
are a common phenomenon. 

One reason why ULXs are so intensely studied is their 
potential connection with intermediate mass black holes 
(IMBHs). Under the assumption that they emit isotropically, 
the Eddington limit gives a lower bound for the mass of the 
central black hole as large as a few hundred solar masses for 
nominal Eddington ratios of ~0.1 — 1. The existence of IMBHs 
would be extremely exciting, as they could be the 'missing link' 
between stellar mass black holes and the supermassive black 



holes in the center of the galaxies (see e.g., lEbisuzaki et alJ 
2001 ). Such objects would have important implications for cos- 
mology and arise in many theories on the collapse of primordial 
stars. 

While there is eviden ce that a few of these objects 
are indeed IMBH s (e.g., IStrohmaver & Mushotzkvl 120031 
ICropper et alJl2004J) . the creation and feeding mechanisms for 
black holes of this mass are still not well understood. The ULX 
phenomenon seems t o be connected t o star formation, see e.g., 
the cartwheel galaxy jGao et al. 2003). This galaxy would need 
hundreds to thousands of IMBHs feeding from a yet unknown 
non-stellar mass reserv oir or the ULX phenomenon is con- 
nected high mas s XRBs jKinj2004 . lGrimm etaD J2002i) and 
iGilfanovl J2004I) show the existence of a universal X-ray lu- 
minosity function that may extent up to ULX luminosities by 
taking the star forming rate into account. There seems to be a 
dependence of the ULX abundance on the galaxy type, since 
dwarf galaxies host many of the nearby ULXs. 

There is also an ongoing discussion whether the inner ac- 
cretion disk temperatures of ULXs are too hot for IMBHs, 
since many have kT,„ ~ 1 keV, which is more similar to thos e 
found in normal XRBs see e.g..lColbert & Mushotzkv ( 1999); 
iMizuno et alJ Jl999l) : iMakishim a et alJ (12000). However, re- 
cent XMM observations find lower inner disk temperatures 
for some ULXs, which is more compatible with IMBHs than 
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earlier resu lts from ASCA observations (see e.g. iMiller et"afl 
l2003[l2004 . 

Besides the explanation that ULXs are accreting IMBHs 
other stellar-mass ULX models have been proposed: 
aniso tropic X-ray emission ("mild bea ming") (|K ing et al. 
20011) , super-Edd ington accretion flows ( Abramowicz et alJ 
1988tlKin gll2004 . and microblazars JKording et alJl2002l) . In 
the microblazar model, the 'ultra-luminous' X-ray emission is 
created by a relativistically boosted X-ray jet. While the jet is 
intrinsically weaker than the accretion disk, Doppler boosting 
can lead to the observed luminosities. If the je ts also emit in the 
radio regime like XRBs in the hard state (cf., Fender 2001) or 
the very high state (see e.g., GRS 1915 Rodri guez et al . 1995), 
the radio emissi on is also boosted and may be observable. The 
jet/disk model jFalcke & Biermam] 1 19951) assumes that the 
accretion flow and the jet form a symbiotic system; if there is 
a disk, there is always a jet in some form. Flat spectrum radio 
emission originates from the self-absorbed part of the jet, while 
one observes a steep synchrotron power law in the X-rays 
(Mai- koff et all 1200 ll) . Thus, one possibility to distinguish 
between the microblazar model and the other possibilities is 
the presence of compact, luminous radio cores at the positions 
of ULXs. 

Galactic XRBs in the very high state are probably the near- 
est cousins to ULXs. Of these GRS 1915+105 can be seen 
as the prototype. GRS 1915+105 shows bright radio flares 
of 1. 5 Jy while the quiesc ent value is approximately 130 
mJy jRodrigu ez et al.l fl995h . An other highly variable XRB 
is Cyg X-3. While this source is also found in the low and 
high state, it shows major radio flares with a flux increase of 
a factor 1 0-100 on a time sc ale of a day which lasts for days 
or weeks (Qgl ev et aljEoOll) . Therefore, it has to be expected 
that also ULXs could be radio transients. These relativistically 
boosted radio flares could be detectable with the VLA. In ad- 
dition, these observations can be used to estimate the number 
of radio flares created in a galaxy. Such estimates are important 
for the design of new digital radio telescopes such as LOFAR 
or the SKA, as they could easily search for transients due to 
their multi-beam capabilities. 

In this paper we describe a systematic search for steady 
state or flaring radio emission from a well defined sample of 
ULXs. In Sect. 2 we present our observation scheme, and in 
Sect. 3 we show our results. Theoretical implications are dis- 
cussed in Sect. 4. Our conclusions are given in Sect. 5. Finally 
we discuss the low luminosity AGN NGC 4736 in the ap- 
pendix. 

2. Observations 

2.1. Sample Selection 

ULXs are defined as having observed 0.3-8.0 keV luminosities 
Lx > 10 39 erg s _1 . We selected nine ULXs from a larger list 
of ~200 ULXs, which were found in Chandra ACIS data that 
had become public by 4 Feb 2003. X-ray luminosities were 
estimated from ACIS count rates using a power-law model with 
photon index F =1.7, and the Galactic absorption column to 
the corresponding galaxy (calculated using the FTOOLS NH 



Table 1. Positions of the ULXs and their 0.3-8 keV X-ray lu- 
minosities. In the fields with only one ULX, the pointing center 
is the ULX. The field of NGC 4736 was centered at the mid- 
dle of both sources. A * in front of the ULX name denotes the 
pointing center. 



FOV 


Name 


RA 


Dec 


X-ray Lum. 
r 39 erg n 


1 * 


M33 


01:33:50.9 


+30:39:38 


l 


2 * 


NGC 2403 


07:36:25.5 


+65:35:40 


2 


* 

3 


M82-D 


09:55:50.2 


+69:40:46 


3 




M82-F 


09:55:51.0 


+69:40:45 


6 


4 * 


NGC 4736-1 


12:50:53.1 


+41:07:13 


1 




NGC 4736-2 


12:50:53.3 


+41:07:14 


1 


5* 


NGC 5204 


13:29:38.6 


+58:25:06 


2 


6* 


NGC 5457-1 


14:03:32.4 


+54:21:03 


2 


7 * 


NGC 5457-2 


14:04:14.3 


+54:26:04 


1 



program). See Tablef^for Lx values for the nine ULXs in our 
sample. 

X-ray point source lists, used to create the larger ULX cat- 
alog, were cre ated using the XASSIS T automatic X-ray anal- 
ysis program llPtak & Griffithsll2003[ URL: www.xassist.org). 
Our nine sample objects are the nearest nine ULXs (D < 5.5 
Mpc) in this list, with DEC > -20°, excluding M81 X-l. The 
latter source was excluded because it has already been stud- 
ied with the VLA (more than 4 hours of integration time, no 
published detections). Nevertheless, they represent a distance- 
limited sample that is unbiased in X-ray luminosity. 

The distance limit was mainly used to reduce the sample 
size. A low distance also helps to give low limits on the to- 
tal emitted radio power of the sources. This selection results in 
nine ULXs in seven fields of views (FOVs): M33, NGC 2403, 
M82 (2 ULXs), NGC 4736 (2 bright point sources, one may be 
the nucleus), NGC 5204, and two fields in NGC 5457 which is 
also known as M101. The positions of the ULXs are shown in 
tabled All host galaxies are spir als except M82, which is cla s- 
sified as dwarf/irregular (see e.g. Jde Vaucouleurs et al!l99ll) . 

2.2. Observing scheme 

The goal of these observations is to search for radio flares and 
continuous emission from ULXs using the VLA. The seven se- 
lected fields were observed eight times each. The observing 
runs were approximately four hours long and occurred between 
June and October 2003. This results in 34 minutes per source 
and epoch including the calibration scans. After the last epoch 
the data of all epochs were combined to give deep observations 
to search for continuous emission. To avoid confusion the ob- 
servations were obtained in the VLA's highest resolution con- 
figurations (A or B configuration). 

The optimal receivers at the VLA to search for weak, flat 
spectrum point sources are the X-band receivers (8.49 GHz). 
They are more sensitive than the C and U band receivers, as 
long as the spectrum is flatter than a spectral index a (S v oc v~ a ) 
of 0.4. We expect that the ULXs have flat or even inverted 
spectrum which could be detectable if they are beamed. Thus, 
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Table 2. Observing Dates of the different Epochs. The sources 
are observed in two groups. The first group consists of M33, 
NGC 2403, and M82, while NGC 4736, NGC 5204, and the 
two FOVs in NGC 5457 belong to the second group. 



Epoch 


Group 1 


Group 2 


Configuration 


1 


2 June 


4 Jun 


A 


2 


3 July 


3 July 


A 


3 


28 July 


5 July 


A 


4 


18 Aug 


17 Aug 


A 


5 


7 Sep 


30 Aug 


A or AnB 


6 


15. Sep 


15. Sep 


AnB 


7 


03 Oct 


07 Oct 


AnB 


8 


16 Oct 


18 Oct 


B 



Table 3. RMS flux values for the different FOVs. Besides the 
RMS of each epoch we give the RMS value for map of the 
combined dataset. All values are given in yuJy. For a discussion 
on M82 see Sect EH 



Epoch 


1 


2 


3 


4 


5 


6 


7 


8 


comb. 


M33 


37 


36 


37 


37 


60 


41 


71 


54 


17 


2403 


31 


30 


32 


31 


59 


38 


45 


37 


13 


N5457I1 


33 


32 


31 


35 


35 


39 


77 


46 


14 


N5457I2 


32 


32 


34 


37 


33 


35 


75 


34 


13 


N5204 


35 


33 


31 


35 


34 


37 


51 


38 


13 


N4736 


45 


47 


49 


55 


45 


61 


63 


62 


21 



we observed at 8.49 GHz using the maximal bandwidth of 
50 MHz. Due to the distribution of the right ascension of our 
sources we had to observe them in two groups (for the groups 
and the observing dates see table|2J. 

As we are searching for very faint radio emission, we 
obtained phase-referenced observations. For all sources VLA 
phase-calibrators could be found within 10°, typically within 
5° of the target source. The cycle time between calibrator and 
source was approximately 7 minutes. The phase correction be- 
tween subsequent calibrator scans was typically < 40° on the 
longest baselines, while the phase corrections on the shorter 
baselines are typically < 10°. Therefore, the phase calibrated 
images should only be degraded by a few percent due to phase 
errors. Amplitude calibration has been done using either 3C286 
or 3C48. The phase-calibrated data has then been imaged with 
natural weighting to achieve the maximal sensitivity. In each 
epoch a source was observed for approximately 21 minutes, 
excluding phase and amplitude calibrator observations. 

For M82, we h ave also reanalyzed archiva l A-configuration 
data presented in lKronberg & Sramekl dl985t) . as they have re- 
ported a radio flare of 7.07 mJy at C-band (6 cm). Intriguingly, 
the position of the flare is near the brightest ULX in this galaxy. 

3. Observational Results 

3.1. Error limits 

The root mean square (rms) flux values of the observed sources 
have been calculated from natural weighted maps. For fields 



with detected sources the largest possible rectangular field ex- 
cluding the source has been used. The theoretical rms value 
for a naturally weighted 21 minute observation with 27 anten- 
nae is, according t o the VLA Observational Status Summary 
dTavlor et al . 2004), 32 //Jy. The rms flux values of the individ- 
ual epochs and fields are shown in table In a single epoch 
the typical observed rms is about 35 fiiy (except for M82, see 
Sect. 13.2. II . If one combines all epochs, the rms goes down 
to * 15/iJy, while the theoretical limit is 11 /iJy. To be able 
to interpret the radio images we have to ask when a peak in 
the map is significant. Chandra has a point spread fu nction 
(PSF) of approximately one arc-second. Wrob eTet alJ (|2001) 
report that astrometric error for phase referenced VLA observa- 
tions (A-configuration) compared to the International Celestial 
Reference Frame as established by VLBI jMa et allll998h is 
less than 10 mas for the X-band. The astrometric uncertainty 
of the VLA is therefore negligible compared to the Chandra 
position errors. 

We know the positions of the ULXs with the accuracy of 
Chandra. During the first five epochs the VLA was in A con- 
figuration and our beam width was * 240 mas. The PSF of 
Chandra which has a full width half maximum of around one 
arcsecond, is therefore covered by 17 beams. The probability 
that we detect a random peak in 17 independent beams is given 
by the error function. Under the assumption that these beams 
are indeed independent, the chance that there is a random 3 <x 
peak at the position of a given ULX is 4 %, the chance of a 4 
cr peak 0. 1 %. We can therefore only accept peaks at the posi- 
tions of the ULXs, if they exceed 4 cr. To search for other radio 
flares in the whole map we have to increase the 4 cr limit. We 
typically map 50" by 50", which corresponds to roughly 14000 
'independent' beams, so there will usually be one 4 cr peak in 
the map. Therefore point sources without known positions are 
not believable unless they have a higher signal-to-noise than 

~5cr. 

3.2. Detections 
3.2.1. M82 

In the starbursting galaxy M82 many compact radio sources 
(CRS s) have been found (see eg., iRronberg & Wilkinson! 
1 19751) . Many of them have been identified as supernova 
remna nts (SNRs). This has been done by direct imaging 
(e.g.. IPedlar et all Il999l) or by spectral observations (see 
lAllen & K ronberg 1998). M82 also shows strong diffuse ra- 
dio emission, making it more difficult to achieve the optimal 
noise limit for the short A-config snap-shot observations (zero- 
spacing, dynamic range limitations). The longer baselines are 
less contaminated by this diffuse emission. Thus, the best rms 
values have been obtained by only using baselines exceeding 
300 k/l for the A and early AnB config observations. This omits 
approximately half of the baselines, increasing the theoretical 
rms by a factor of =s V2 to 44 /jjy. The rms values for a single 
epoch still do not reach this value but are a bit higher: « 70 yuJy. 
For the 6th to the 8th epoch smaller cut-offs had to be utilized 
while the array configuration approached B-configuration. In 
table|4]we give the X-ray and radio positions. While the X-ray 
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PLot file version 1 created 24-NOV-2004 13:20:47 
CONT:M82IXOS IPOL 8460.100 MHZ M82-1-6.ICL001.1 




09 55 51.5 



51.0 50.5 50.0 

RIGHT ASCENSION (J2000) 
Cont peak flux = 1.7657E-02 JY/BEAM 
Levs = 6.000E-04 * (-10, -9, -8, -7, -6, -5, -4, 
-3, -2, -1,1,2, 3, 4, 5, 6, 7, 8, 9,10) 



49.5 



Fig. 1. Map of M82 around the ULXs D and F (see table|4}. The brightest radio source in the field is e->_ which is coincident with 
X-ray point source E. The FWHM of all X-ray point sources have been set to 1". Even though the radio source d\ is outside of 
the r=0.5" circle for ULX D, its offset (0.77") is consistent with the ULX X-ray position. 



positions are named with capital letters, the corresponding ra- 
dio detections are named as indexed small l etters. The X-ray 
sources D and F are the well known ULXs ( Matsumoto et al. 




Our VLA observation was centered on ULX D of table 0] 
Thus, the bright X-ray point source A is 18.8" away from the 
phase tracking center. As we are observing in continuum mode, 
our sensitivity is already reduced for this ULX due to band- 
width smearing. The distance corresponds to 72 beam widths , 
thus, the peak flux will be reduced by 5% JTavlor et al . 2004). 
We have not detected any radio emission from source A. 

The brightest CRS in the FOV (e 2 , 41.95+57.5 in the no- 
tation of Mux low et alj fl994. abbreviated to MPW) is in good 
agreement with the position of a brigh t X-ray point source (E) . 
This source is probably a SNR (e.g., Allen & Kronberg 1998) 
that may be immersed in a molecular cloud ( Ped lar etalJl 999). 
As the radio and X-ray positions coincide, the radio and X-ray 
maps are correctly aligned. 

Near the ULX D we have detected a CRS (di) with a contin- 
uum flux of 1.8 mJy in all eight epochs observed in this cam- 
paign. This source is also detected by MP W with 1.8 mJy a s 
41 .3+59.6 and is identified as a SNR jAllen & Kronberdl9 98). 
As expected, the flux of the source is stable during the 8 epochs 
within lcr of the rms. It has a distance of 0.77" to the ULX, 
well in agreement with the Chand ra position. We did not detect 
the flaring source of Kronberg & Sramekl Jl985t) in our VLA 
observations. 



For the ULX F we find two possible radio counterparts as 
shown in tabled The brighter one of both is also given in ear- 
lier papers as 42.2 1 +59.0 while the weaker one is barely visible 
in the maps of MPW but not given in their table. Both sources 
are not variable over our observing period. 



Besides the CRS near the reported positions of the 
ULXs we have not f ound any radio flare in M82. However, 
iRronberg & Sramekl {l985) found a 7.07 mJy flare near the 
brightest ULX (D) in M82 (2.3 x 10 19 W/Hz). We have re- 
analyzed their data for direct comparison with our results. The 
position of the M82-flare is is 0.56" away from the Chandra 
position of the ULX, see table |4] It is already visible in the 
phase-referenced images before any self-calibration. In our re- 
analyzed maps the flux of the fla ring source is slightly highe r 
(» 9 mJy) than the value found by Kronberg & Sramek (1985). 
The brightest CRS at that tim e (e?) had a flux of 105 m Jy in 
1981 (108 mJy as measured bv lKronberg & Sr amek 1985J), cur- 
rently this CRS has an X-band flux of 12.3 mJy. The M82-flare 
is stable throughout the 7 hour observation in 1981, and is also 
visible if one divides the observation into 1 hour blocks. Even 
though, the detected radio flare is near the ULX position, it 
can well be, that this source is also a SNR. The radio decay 
timescale and the spec trum is similar to the SN 1983n in M83 
jKronberg et all2000l) . which was identified with a SN associ- 
ated with a massive star ( Sra mek et alJl 984). 
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Table 4. Radio detections near detected X-ray point sources in 
M82. A flux of 1 mJy corresponds to 3 x 10 18 W/Hz at the 
distance of M82. 



Name 


X-ray Position 


Radio F. 


Dist. 




Radio Position 


LmJyJ 


arc sec 


A 


09:55:46.59 +69:40:40.9 






B 


f\C\ CC A £L T 1 , f l\ A f\ T7 C\ 

09:55:46.71 +69:40:37.9 






C 


f\f\ C C \ ^1 AO / t \ A 1 \ C t \ / 

09:55:47.48 +69:40:59.6 






D 


09:55:50.15 +69:40:46.5 






d\ 


f\C\ CC C f\ f\C . /Tl A f\ AC f\ 

09:55:50.05 +69:40:45.9 


1.8 


0.77 


M82-flare 


r\c\ c c c r\ 1 r\ ,/"n a f\ a f f\ 

09:55:50.19 +69:40:46.0 


9.0 


0.56 


E 


r\c\ cc c t\ q/i , / i\ a r\ ho o 

09:55:50.70 +69:40:43.8 






e\ 


09:55:50.39 +69:40:44.3 


1.5 


1.7 


(-'2 


09:55:50.69 +69:40:43.8 


12.3 


0.06 


ei 


09:55:50.82 +69:40:44.7 


1.6 


1.08 


F 


f\C\ Cc c 1 m , S~ (~\ A f\ AC 1 

09:55:51.03 +69:40:45.1 






^3 


09:55:50.82 +69:40:44.7 


1.6 


1.21 


A 


09:55:50.95 +69:40:45.3 


1.9 


0.47 


G 


09:55:51.29 +69:40:43.9 






gi 


09:55:51.22 +69:40:44.5 


1.3 


0.67 


g2 


09:55:51.29 +69:40:44.3 


1.3 


0.37 


83 


09:55:51.41 +69:40:44.5 


0.4 


0.87 


H 


09:55:51.48 +69:40:36.0 






I 


09:55:53.43 +69:41:02.0 
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Fig. 2. The double source in NGC 4736. The upper source is 
called source a and the lower and stronger source is called b. 



Table 5. Other radio detections 



Name 


Radio Position 


Radio Flux [mJy] 




NGC 4736 




N4736-a 


12 50 53.03 +41 07 13.6 


1.6 


N4736-b 


12 50 53.08 +41 07 13.0 


1.2 




M101 




MIOl-a 


14 03 31.31 +54 21 14.9 


0.092 


M101-6 


14 04 14.32 +54 26 09.4 


0.057 


MIOl-c 


14 04 14.25 +54 26 09.4 


0.060 


Low significance flux peaks, see discussion 


MlOl-o 


14 04 14.21 +54 26 02.6 


0.15 


N2403-Q- 


07 36 25.49 +65 35 39.8 


0.14 



3.2.2. NGC 4736 

During all epochs we detect a double radio source near ULX 
NGC 4736-1 as shown in Fig. [2 Their separation is about 1", 
which corresponds to 20 pc at the galaxy for a distance of 4.3 
Mpc. The position of the stronger radi o source coincides with 
the position of the LLAGN reported bv lNaear et alJ J2002I) and 
is in agreement with the Chandra position of ULX NGC 4736- 
1 (see tableland Fig . - No radio emission is detected at the 
position of the other ULX NGC 4736-2 (see Fig.0. As one of 
the two X-ray sources is probably the nucleus, we assume that 
the stronger radio source of the two (source b) is the nucleus, 
while the other radio source is of unknown nature. As this dou- 
ble structure is probably not related with the ULX phenomenon 
we discuss the findings in the appendix A. 



3.2.3. NGC 5457 (M101) 

In the combined maps of both fields of view the galaxy 
NGC 5457, which is also known as M101, in we find emission 
at a 5 <x level (Ml 01 -a), which is not at the positions of the 
ULXs. At the same position is a 4cr peak in the eighth epoch. 
If one excludes this epoch from the map the source is still visi- 
ble but has a reduced flux of 72.9 //Jy. In the second FOV in 
NGC 5457 there seems to be a double source: MlOl-b and 
MIOl-c. Both sources are 5cr detections. The nature of these 
sources is unknown. 

3.3. Non-detections 

There are no significant detections found in M33, NGC 2403, 
and NGC 5204 in the individual epochs or the full data with a 
signal to noise higher than 4. In M 1 1 we do not find significant 
detection in the individual epochs, for continuous emission see 
Sect. 13.2.31 Here we limited the search radius to 2" from the 
ULX Chandra X-ray positions. 

In NGC 2403 we found in the fourth epoch a peak in the 
map (N2403-Q') which is a 4.4 sigma detection. The VLA po- 
sition coincides with the ULX position within the accuracy of 
Chandra. However, the emission seems to originate only from 
a short time period (5 of 21 minutes). If these times are omit- 
ted, the peak in the map vanishes. Some XRBs show such short 
time radio variability (Ogley et al. 2001). No second short radio 
flare was seen in our data. So, it is not possible to confirm such 
a short flare, especially as the significance is only marginally. 
We therefore do not consider this 'flux peak' as a detection of 
a flare. 
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Table 6. The radio power of our sample for continuous emis- 
sion and flares. The radio detections in M82 are probably SNRs 
while in NGC 4736 we probably detect the nucleus of the 
galaxy. For all other sources we present upper limits for con- 
tinuous e mission over a ll epochs and flares. The distances are 
taken from lTullvl(ll988l) . 



Name 


Dist. 




L R,flaie 


^R,cont 




[Mpc] 


[io- -Ij 


[1033S] 




M33 


0.7 


1 


< 0.43 


< 0.2 


NGC 2403 


4.2 


2 


< 13 


< 5.4 


M82-D 


5.2 


3 


<45 


290 


M82-F 


5.2 


6 


<45 


307 


NGC 4736-1 


4.3 


1 


< 21 


177 


NGC 4736-2 


4.3 


1 


< 21 


<9.2 


NGC 5204 


4.8 


2 


< 18 


< 7.1 


NGC 5457-1 


5.4 


2 


< 23 


<9.1 


NGC 5457-2 


5.4 


1 


< 23 


<9.1 



In M101 we detect a 4.2 cr peak (MIOl-a, see table in 
the second VLA pointing centered on the ULX NGC 5457-2. 
The radio source is 1.4" away from the Chandra position of 
the ULX, further away than the typical position accuracy of 
Chandra, but it is still possible that it is connected to the ULX. 
However, besides this peak there is another similar maximum 
another arcsecond away, and both maxima are found on a side- 
lobe. Thus, it is possible that the source could be an instrumen- 
tal artifact. Besides this, we do not detect significant (> 4<x) 
radio emission from the ULXs in the individual epochs. 



4. Theoretical Interpretations 

4.1. Flares 

The nearest known Galactic cousins of ULXs are the highly 
accreting black hole XRBs. XRBs in the low hard or the high 
state are highly va riable in the radio and X-ray regime (c.f., 
Ivan der Klisl Il989). The radio flux can vary by a few hun- 
dred percent, for example Cyg X-l has a typical 15 GHz ra- 
dio flux of approx imately 10 mJy and reaches up to 30 mJy 
JPottschmidt et aljl2003l) . As already mentioned in the intro- 
duction stronger accreting objects like GRS 1915+105 or Cyg 
X-3 show violent radio flares of a factor 10 to 100. During 
these flares GRS 1915+105 reaches 1.5 Jy IRodriguez et alj 



1995) and Cyg X-3 goes up to more than 15 Jy (Watanabe et al. 
1994), while the 'steady' state emission is around 100 mJy for 
both sources. The brighter radio bursts of Cyg X-3 may be due 
to a higher Doppl er factor, as the inc lination of this source is 
only 14° jMioduszews ki et alJl200ll) . GRS 1915+105 on the 
othe r hand is seen nearly edge on (angle to the line of sight: 
66° . IRodriguez et alJl995l) . 

One possibility to estimate the luminosities of radio flares 
one could expect from beamed emission of ULXs is by tak- 
ing the transients GRS 1915+105 and Cyg X-3 as an example. 
The average distance of our ULXs is approxim ately 4.6 Mpc . 
GRS 1915+105 has a distance o f 1 1 kpc JFender et aljfl9 99b) 
while Cyg X-3 is 9 kpc away dPredehl et alJl2000l) . Without 



any relativistic beaming these sources would have a flux of 8.5 
yuJy and 57 //Jy, which is below our detection limit. However, 
already a mild beaming factor of 20 for GRS 1915 or an addi- 
tional factor of 3 for Cyg X-3 would bring the flares into our 
detection limit (0.15 mJy). As GRS 1915 is seen with an in- 
clination angle of 66°, its Doppler factor will be around one. 
Thus, a beaming factor of 20 can be reached with a moderate 
Lorentz factor of T — 3 ( Lind & Blandfor dl l!985l) . if the jet 
points roughly at the observer (inclination angle < 15°). The 
Doppler factor of Cyg X-3 is uncertain, as the observed lumi- 
nosity is probably already beamed. But an additional factor of 
3 is easy to obtain if the jet is pointing directly to the observer. 
If the ULXs are flaring in radio we should be able to detect the 
beamed radio flares. 

We have not found a single significant flare for all our 
sources in the monitoring campaign. In the individual epochs 
our 4<x sensitivity is on the average 0. 15 mJy. 

T he distance of M33 is around 0.84 Mpc ( Freedm an et alJ 
119911) . Besides this nearby galaxy all other observed galaxies 
have distances in the range from 3.6 Mpc to 5.4 Mpc. Thus, the 
detection limits of flares will only vary by a factor of two for 
those galaxies. For the average distance of 4.6 Mpc the upper 
limit on the radio power S v of flares is 3.8 x 10 17 W/Hz. This 
corresponds to a 5 GHz radio luminosity (vS v ) of 1.9 x 10 34 
erg/sec assuming a flat spectrum. For M33, however, the limit 
is reduced to a radio power of 1.1 x 10 16 W/Hz and a radio 
luminosity of 5.7 x 10 32 erg/sec. 

The non-detection may be due to the unknown time scales 
of the radio flares in ULXs. We observed our sources once or 
twice a month. If the time scale of a radio flare is only a day, 
we are strongly under-sampling the radio light-curve. The time 
scale of the boosted flares is unclear and depends on the physi- 
cal process creating the flare. If the flare is created inside the jet, 
e.g., similar to the shock in jet models used forblazars (see e.g., 
Marsche r & Gearll 1985b , the observed time scale of the flare 
will be reduced by the Lorentz factor. For F x 5 the time scale 
could be as short as a few hours. However, it will be extremely 
bright, even a few mJys are possible. On the other hand, if the 
flare is created by enhanced injection of material into the jet by 
the disk for an extended time, the observed time scale will be 
the same as the intrinsic one. 

The fact that we have not found a single flare in all our 
epochs can be translated to an upper limit for the probability 
that an average ULXs of our sample has a radio flare brighter 
than 3.8 xlO 17 W/Hz, which should be detectable. If we assume 
that the probability that an ULXs is flaring at a given time is 
similar for all observed ULXs, and one flare is uncorrelated to 
earlier or later flares, the flares should be Poisson distributed. 
The probability distribution describing how many events we 
detect given the probability of a detection is 

exp-' 1 A" 



P(A, n) 

n: 

where A denotes the expectation value of the distribution and n 
is the number of events. 

For this study we have to exclude M82 as the rms in the 
M82 maps is much higher than in the other fields. The other 
six fields have been observed for eight times which yields 48 
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samples. Let p denote the upper bound of the duty cycle of 
the ULX. Thus, it is an upper limit for the probability that one 
ULX flares at a given moment, p has to be chosen, such that we 
should have almost certainly detected at least one event if the 
probability that a given ULX flares is p. Here p will be chosen 
such that we should detect one or more flares with a probability 
of 95%. This leads to 

n=co 

0.95 = Yj p ( 48 P' ») = 1 - exp~ 48p . 

This results in p — 0.06, i.e. the duty cycle of radio flares ex- 
ceeding our detection threshold in ULXs is < 6%. The expec- 
tation value for the number of flares is 2.9, our detection of no 
flares is in agreement with this value at a 5 % level. 

If we know the time scale of a typical radio flare in an ULX 
we can convert the upper limit of the flaring probability to an 
upper limit of the number of flares a ULX can have in a year. 
Let tffpiare denote an average length of a flare. As the radio flares 
of GRS 1915 have a time scale of days, we will use <5fFiare ~ 2 
days as a reference. Thus, the upper limit of the probability to 
detect a flare in a single observation corresponds to an upper 
limit of 1 1 radio flares per year. This upper limit is still 
higher than the number of bright flares in GRS 1915+105. 

Besides looking for radio emission from ULXs we can use 
these observations to derive upper limits on the amount of ra- 
dio transients happening in a galaxy. Our galaxies have an av- 
erage distance around 4.6 Mpc and we mapped an area of 51", 
this results in an observed area of 1.3 kpc 2 . We are searching 
for radio flares from unknown sources, therefore, it is also not 
known which is the quantity correlating with the number of 
flares. While for ULXs the star formation rate might be a good 
quantity it could simply be the number of stars or black holes 
for an other class of flaring objects. As we are observing a cou- 
ple of small fields in different galaxies, it is very hard to derive 
the exact number of observed solar masses or even the observed 
star formation rate. As most of our observed galaxies (besides 
M82, which we exclude) are non starbursting galaxies, the star 
formation rate per solar mass will be of the same order of mag- 
nitude. In order to get a rough estimate how much mass we have 
actually observed, we assume a similar mass density as in our 
Galactic neighborhood. In our Galactic neighborhood the mass 
density of the disk is of the order of 200 M Q per pc 2 . Therefore 
one field of view observes a mass of approximately 3 x 10 8 M . 
Thus, we expect less than 0.2 flares in a single observation of 
10 9 Mo. 

4.2. Background sources 

We have detected continuous radio sources only inside 
NGC 5457 (excluding M82). In all other fields there were no 
backgroun d sources visible. How ever, we only map an area of 
51" x 51". lFomalont et all (J2002) find that the number density 
of sources above a flux density S is approximately 

/ „ >-l.ll±0.13 

N = 0.099 ± 0.01 ( arcmhr 2 . (1) 

\40/Jy/ 

Our limiting sensitivity (4cr) is approximately 56 /iJy. Thus, 
this predicts that there should on average be 0.049 ± 0.007 



sources in a given observed 51" x 51" map. Our findings are in 
agreement with this result. The probability that a background 
source is found coincident with a ULX, within the position ac- 
curacy of Chandra (1 arcsec), is below 2 x 10~ 5 . 

4.3. Steady State Emission 

In M82 two reported Chandra ULXs positions are in agreement 
with the positions of CRSs, that are probably radio SNR. For all 
other galaxies we have not detected any continuous emission 
near the ULXs. The upper flux limit for the observed sources 
is * 60yt/Jy (4cr). For the steady state emission ULXs should 
lie on the Rad io/X-ray correlation, which was first established 
for XRBs JCorbel et alJ l2000l 120031 iGallo et all l2003t) and 
has b een extended to AGN jMerloni et alJl2003l iFalcke etail 
|2004 . Simple jet scaling (see e.g.. lFalcke et alJl2004 predicts 
that the X-ray and Radio emission are correlated and scale ac- 
cording to: 

L x cc 4 38 M - 81 (2) 

where Lx and Lr denote the X-ray and radio luminosities. 
The exact value of the exponents depends on the spectral in- 
dices of the Radio and X-ray emission. In Fig. [5] we plot the 
3-9 keV X-ray luminosity against the 5 GHz radio luminos- 
ity. The radio luminosity is always calculated assuming a flat 
spectrum, i.e., it is defined as L$ qjj z = 5GHzS 5 qjj z where 
$5 GHz denotes the radio power at 5 GHz. The Chandra X- 
ray luminosities are corrected for the different energy band as- 
suming a photon index of 1.6. Besides a sample of low hard 
state XRB (GX339-4, V404 Cyg, XTE Jl 1 18+480, 4U 1543- 
47, GS 1354-64) we show the two XRBs used to extrapolate to 
the expected radio flare luminosities: GRS1915 an d Cyg X-3. 
The data-points have been taken from Gall o et alJ (12003). The 
brightest 1.5 Jy flares of GRS 19 15 are not included in the plot. 
Additionally we present a sample of AG N consisting of low 
luminosity AGN, FR-I radio galaxies ( Fan aroff & Rile\ll974l) . 
and Radio and X-ray selected BL Lac objects. This figure is 
similar to Fig. 3 in Falc ke et alJ 12004). Furthermore, we show 
the theoretical radio/X-ray scaling for Black holes of different 
black hole masses. 

A strict radio/X-ray correlation will only be valid for 
the steady state emission of low/hard state, jet-dominated 
sources. For high state objects the radio emission wi ll be 
quenched (see, e.g. JFender et al.ll999allGallo et alJ2 003). For 
one possible cor relation extending to the very high state see 
Maccarone et alJ ( I2003I) . Furtheremore, we can not expect that 
the correlation holds for the flares of the transient sources. 
Those flaring sources are only included for comparison. 

The empirical correlation can be used to compare the upper 
limits (4<x) of this radio monitoring campaign with the expec- 
tations. The radio flux of high or very high state objects may be 
reduced compared to the low/hard state objects. The correla- 
tion gives an upper bound for the expected radio flux. Only the 
upper limit for M33 is close to the value expected from the cor- 
relation for stellar mass black holes. All other upper limit are 
far above the expected radio flux from microblazars with stellar 
mass black holes. Our non-detections are consistent with either 
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27.5 30 32.5 35 37.5 40 42.5 

Lg L R / (erg/s) 

Fig. 3. The upper limits of the ULXs in respect to the low/hard state XRB and AGN population. The lines correspond to different 
black hole masses. From Top to bottom 6, 6 x 10 2 , 6 x 10 6 , 6 x 10 9 M o . The left ULX radio detection represents NGC 5408, while 
the right is the radio flare in M82. The correlation should not hold for fares, this was just plotted for comparison. As GRS 1915 
and Cyg X-3 are both radio transients, they are not expected to follow the correlation. 



beamed or non-beamed radio emission from stellar or interme- 
diate mass black holes that are < 1O 3 M . 

4.4. Other radio detections of ULXs 

Up to now there are two radio detections mentioned in the lit- 
erature, here we discuss their int erpretation within the microb- 
lazar model. iKaaret et ail (Eo 03) have detected radio emission 
from an ULX with an X-ray luminosity 1.1 x 10 40 erg/sec in 
NGC 5408. They find a 4.8-GHz radio flux of 0.26 ± 0.04 mJy 
well above our detection limit (17 cr for continuous emiss ion). 
Assuming a distance of 4.8 Mpc jKarachentsev et alT2 002) this 
corresponds to a radio power 7.2 x 10 17 W/Hz, which is also 
above the upper limits for our sample of ULXs (1.5 x 10 17 
W/Hz). For the X-band they give an upper limit of 0.12 mJy 
(3 cr). The spectral index is therefore larger than a > 1.0. We 
would expect that the radio cores of microblazars have a flat to 
inverted spectrum, by far flatter than the observed steep spec- 
trum. The radio source does not seem to be a transient as it 
was detected with a roughly similar flux in a subsequent obser- 
vation (Kaaret priv. comm.). As we have the radio and X-ray 
luminosity we can use the radio/X-ray correlation to compare 
the ULX in NGC 5408 with XRBs and AGN. The source seems 
to be too radio loud for a stellar mass black hole, as shown in 
Fig-El Under the assumption that the source should lie on the 



correlation, we can derive a mass for the object of ~ 1000M o . 
Otherwise, the observed emission could be explained by a flare, 
or by extended emission like a radio supernova. At least this 
object seems to be different than the ULX sample observed by 

us. 

iKronberg & Sramekl \ 19851) found a 7.07 mJy flare in M82 
(Lr = 1.1 x 10 36 erg/sec ), which we also confirmed by re- 
analyzing their data. The flare position is within the errors of 
the brightest ULX in M82 which has up to 9 x 10 40 erg/sec. 
This flare is bright compared with GRS 1915, which we ex- 
pect to show boosted 0.4 mJy flares. However, the value does 
not seem to be unrealistic for a major flare and the uncertainty 
of the flux estimate. Up to now no second event has been re- 
ported, even though M82 has been frequently observed. Thus, 
these events seem to be extremely rare. This is consistent with 
our conclusion that the duty circle of ULX flares in our sam- 
ple is very low (< 6%). For comparison we have also shown 
the flare in Fig. [5] However, the radio/X-ray correlation is only 
valid for steady state emission. 

5. Conclusions 

We have monitored nine ULXs during five months with the 
VLA. Our limiting sensitivity (4 cr) of the individual epochs is 
« 0.15 mJy (4x 10 17 W/Hz for average distance, rms « 36yuJy), 



Elmar K6rding,Edward Colbert and Heino Falcke: A Radio Survey of ULXs 



9 



with the exception of M82. We have not found any significant 
radio flare in the sample. If they last only a few days, we are 
heavily under-sampling the light-curve. However, we can give 
an upper bound of 6% for the probability to find a flare brighter 
than 4 x 10 17 W/Hz in a given observation of an ULXs. This 
translates to an upper bound of 11 ^ flares per year with a 
time scale dt^e- This is much more than observed for example 
in GRS 1915. Thus this monitoring campaign can not rule out 
or strengthen the idea of relativistically beamed flares in ULXs. 

We have reconfirmed the 7 mJy flare near the brightes t 
ULX in M82 in archival data from Kronbe r g & Sram ek (1985). 
This corresponds to a radio power of 2.3 x 10 19 W/Hz. Such a 
radio flare can be well understood in the context of relativisti- 
cally beamed emission from the UL X. An o ther possible expla- 
nation is a SN similar to SN 1983n (Kronbe rg et alHo OO). 

We have detected a yet unreported variable double source 
in NGC 4736. As the positions of the stronger source coincides 
with the position o f the nucleus of the LLAGN as reported in 
Nag ar et all (Eo02). we assume that this double source is prob- 
ably connected to the LLAGN and not to the phenomenon of 
ULXs. 

The search for continuous emission from ULXs was only 
successful for M82 (limiting flux: as 60 /iJy, Act). Here two 
CRSs (probably SNRs) could be associated with ULXs. Thus, 
ULXs do not generally emit steady-state radio emission above 
radio powers of 1.5 x 10 17 W/Hz. Using the Radio/X-ray cor- 
relatio we have shown that our non-detections are consistent 
with either beamed or non-beamed radio emission from stel- 
lar or intermediate mass black holes that are < 1O 3 M . The 
search for a radio loud ULXs is therefore still open. Our re- 
sults suggest that detections will be extremely rare for current 
telescopes. Nevertheless, the investigations of these intriguing 
sources have to continue. 
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Appendix A: The LLAGN in NGC 4736 



As we have already mentioned in Sect. 13.2.21 we have detected 
a double source in NGC 4736 show in Fig. |2] One of these 
sources is stronger (2.2-1.2 mJy, source b) and variable while 
the other source seems slightly extended and is fairly stable 
in flux (1 mJy, source a). The positions are given in table [5] 
The peak fluxes and the ratio of the two sources is shown in 
table lA.ll for a light-curve of the peak fluxes see Fig. IA.ll 

Due to the low fluxes of both sources it is not possible 
to self-calibrate the data itself. The fluxes can therefore be 
reduced by phase errors. They can furthermore be changed 
by amplitude calibration. However, the flux-ratio of the two 
sources should be independent of these effects. In the radio the 
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Fig. A.l. Light-curve of NGC 4736. Note that the VLA array 
configuration changed from A to B during the last epochs. 



variability can also be seen, it varies between 2 and 1.2 on a 
time scale of months as also shown in the table. We therefore 
conclude, that at least the stronger source b is variable. This 
indicates that source b is a compact object, probably the nu- 
cleus of the low-luminosity AG N. It h as the position that was 
previously published by iNagar et aU (|2002l). Source a i s not 
mentioned in earlier publications, but Na gar et alJ (|2002) rms 
value was significantly higher (0.34 mJy). We have reanalyzed 
their data, and detect a 4 cr peak at the position of source a ( 1 .4 
mJy). Source a was therefore already visible in March 1998. 

For the discussion of the other source we have to consider 
that the array configuration of the VLA changed between the 
epochs five and eight from A-configuration to B-configuration, 
increasing the beam size from 260 mas to 840 mas. Extended 
emission will therefore contribute more to later epochs than to 
the first. This explains the increasing peak fluxes of the source 
a from 1 mJy to 1.7 mJy in the last epochs. The effect that the 
peak flux increases with the beam-size can be seen if one uses 
a uvtaper of 400 kA, which roughly doubles the beam-diameter. 
The peak flux increases from 1.1 to 1.3 mJy in epoch one. The 
extension of source a is also visible in the difference between 
the peak and the integrated flux. For source b these two values 
are usually similar, while for source a the integrated value is 
higher by approximately 30 %. 

It is puzzling that the integrated fluxes are more variable 
(1.1 - 1.8 mJy) than the peak fluxes. We checked that this effect 
is not due to phase errors, epoch one and two have phase jumps 
between successive calibrator scans of less then 40°, usually 
less than 10°. This would result in only a few (< 6%) ampli- 
tude loss. One problem of the integrated flux is that the inte- 
gration boxes are small, as both point sources are less than an 
arc-second away. Another possibility could be that source a is 
barely resolved, small changes in the beam size and the atmo- 
sphere could than lead to the observed variability. An extended 
(> 200 mas) source which is variable on time scales of months 
at a distance of 4 Mpc is not possible, as its size would be larger 
than 5 pc. However, there could be a compact substructure in 
the second source which is varying. 

During the last epoch we also obtained a C-band map of 
NGC 4736. As the VLA was already in B-configuration the two 
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sources could not be separated. Both sources combined give a 
peak flux of 4. 1 mJy and an integrated flux of 7.3 mJy. 

To check the spectrum of the two sources we have also ob- 
tained U-band images during the last epoch. Source b has a 
peak and integrated flux flux of 1.7 mJy. The weaker source 
two has a peak flux of 1.1 mJy and 0.95 mJy integrated flux. 
The missing integrated flux is probably due to the small in- 
tegration box around the peak needed to avoid source b. The 
diameter of the beam was 480 mas. We have reanalyzed the old 
data from Mar 1998 to check for the second source. Both peaks 
are visible in the map with fluxes of 2.235 mJy and 1.403 mJy 
respectively. However, the rms was 0.34 mJy making accurate 
flux measurements impossible. 

We can not compare the U-band flux directly to the X- 
band flux, as the beam size of the X-band observations is about 
twice the U-band beam size. As the flux of the second source 
is roughly constant during the first epochs, we compare the U- 
band flux with the last epoch in A-configuration and taper the 
data to get a beam-size of 480 mas. The flux of the second 
source increases to 1.41 mJy. This results in a spectral index 
between 15 and 8.4 GHz of a = 0.42. 

To check whether these sources are compact we have ob- 
served NGC 4736 with the European VLBI Network. While 
source b was detected with a correlated flux of 0.86 mJy source 
a was not detected. The rms in this observation was 0. 1 mJy. 

No other low- luminosity AGN (LLAGN) in 2 cm survey of 
iNaear et alJ J2002I) shows such a double structure. The nature 
of the second source is unknown. Ma oz et alJ (11995) report for 
this gala xy a second UV s ource 2.5" away from the nucleus 
(see also Maoz et al. 1200.4 and suggest that this galaxy could 
be in the final stage of a merger. We have not detected their 
second source. While the position angle of the UV source is 
-3° the second radio source is at -49°. We conclude that source 
a could either be a hotspot or a jet feature of the LLAGN. The 
spectral index is fairly flat. Thus, the hotspot has to be self- 
absorbed. If this source would be observed at a greater distance 
and the two components unresolved, the spectrum would peak 
at a few GHz, as source b has a flat spectrum. This source could 
therefore be a close-up of a GHz peak source (GPS). An other 
explanation is that inside the extended emission there is a com- 
pact flat spectrum source below the detection limit (0.3 mJy) of 
our EVN observations. If confirmed, this would be spectacular 
as this could either be a second nucleus or the radio core of an 
ultra-luminous X-ray source. 
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